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Abstract 

The resummed transverse momentum distribution of the Higgs boson in gluon fusion 
through LO+NLL for small transverse momenta is considered, where the Higgs is 
produced through a top- and bottom-quark loop. We study the mass effects with 
respect to the infinite top-mass approach. The top-mass effects are small and the 
heavy-top limit is valid to better than 4% as long as the Higgs' transverse momentum 
stays below 150 GeV. When the bottom loop is considered as well, the discrepancy 
reaches up to about 10%. We conclude that bottom-mass effects cannot be included 
in a reasonable manner by a naive reweighting procedure in the heavy-top limit. We 
compare our results to an earlier, alternative approach based on POWHEG. 

1 Introduction 

The recent discovery of a scalar particle [l|[2], which may turn out to be the standard 
model (SM) Higgs boson, is one of the biggest achievement of particle physics in the 
last years. The next task is to verify that this particle brakes the electroweak symmetry 
in a certain theory, e.g. the SM, so that massive particles obtain their masses through 
the Higgs mechanism. For this purpose its fermionic and bosonic couplings need to be 
determined, using precision predictions from the theoretical side fs] and comparing them 
to experimental data. Besides the total cross section it might be helpful to use differential 
quantities to disentangle the Higgs couplings to the various particles. 

The most important production mechanism of the Higgs boson in the SM is the gluon 
fusion process, where the Higgs-gluon coupling is mediated through a top-quark loop. 
Higher order corrections are usually calculated in the so-called heavy-top limit, which 
is an effective theory approach, where the top quark is assumed to be infinitely heavy. 
Calculations that keep track of the full top-mass dependence can usually be performed 



only at one perturbative order lower than in the heavy-top approximation. The uncertainty 
constituted by this approach is very specific to hadronic Higgs production in the SM. It has 
been shown for the total inclusive cross section at next-to-next-to-leading order (NNLO) 
[4j|6] that the heavy-top approximation is valid to better than 1% [7-10 . Rather few 
studies aim to validate the heavy-top limit for distributions 11-14 . Generally speaking, 



it was found that the heavy-top limit works well as long as the transverse momentum of 
the Higgs is below the top mass. 

Another specific uncertainty of hadronic Higgs production in the SM emerges from the 
bottom- loop contributior[^ Although suppressed by their couplings, in gluon fusion the 
Higgs-gluon coupling can be mediated by any quark loop. While the effect on the four 
lightest quarks q S {u, d, s, c} is negligible, the bottom quark has a sizable contribution to 
hadronic Higgs production. There is no effective theory approach for the bottom quark 
feasible with current technology. Thus for the bottom-quark contribution one has to stick 
to the perturbative order where the bottom loop can be included. Besides that, one 
approach which might be followed is to account not only for top- but also for bottom-mass 
effects by including higher order corrections through reweighting of the mass effects at 
lower order by the heavy-top limit. We will argue that for certain quantities this appears 
not to be a good approximation. 

Similar to the top-quark effects, rather few studies quantify the importance of the 6-loop 
contribution in the SM 113-15 and explore the best way to include it [16-18 . In the case 



of the total cross section, for example, it contributes 7% at next-to-leading order (NLO) 
for a Higgs mass of 125 GeV at the Large Hadron Collider (LHC) with 8 TeV collider 
energy [15yi9| and cannot be approximated by the heavy-top limit. A similar study to 
the one in this paper has been done in Ref . [l4j , where the matched parton shower in the 
POWHEG framework |20| with full top- and bottom- mass dependence is studied. 

One of the most important differential observables is the transverse momentum distribu- 
tion of the Higgs boson. In the effective theory approach of the gluon fusion process it 
has been calculated numerically [21] and analytically [22] at NLO. Furthermore, the fully 

These calculations are 



differential Higgs cross section is known through NNLO [23 , 24 



only valid for sufficiently high transverse momenta, because they diverge logarithmically 
in the limit where the transverse momentum of the Higgs vanishes. To obtain a reliable 
prediction in the small transverse momentum region, these logarithmically enhanced terms 
need to be resummed to all orders in perturbation theory. This has been worked out in the 



heavy-top limit at leading logarithmic (LL), next-to- leading logarithmic (NLL) [25 -27 and 



next-to-next-to- leading logarithmic accuracy (NNLL) [28 - 30 in the case of the gluon fusion 
process. The matched NLO+NNLL transverse momentum distribution can be calculated 



with the publicly available program HqT 28 -30 



In this paper we show results for the transverse momentum spectrum of the Higgs boson 



^Throughout this paper we consider the interference terms of the top- and bottom-quark amphtudes 
as part of the bottom-quark contribution. 



2 



in the gluon fusion process with exact top- and bottom-mass dependence, including the 
resummation of large logarithmic terms for small transverse momenta at LO-I-NLL. We 
compare it to the effective theory approach and study the quality of the heavy-top limit 
to describe the exact top-mass cross section. Furthermore, we quantify the impact of 
the bottom-loop contribution in the SM. We argue that the bottom-mass effects should 
be included at the order where their calculation is feasible and not by using a naive 
reweighting procedure in the heavy-top limit. As a result, we find that the effective 
theory approach approximates the exact top-mass dependence at the level of 0.5% in the 
region where resummation is important and to better than 4.5% for transverse momenta 
below 150 GeV. The uncertainty induced by a missing bottom-quark contribution is of the 
order of 10%. Finally, we compare our results for the analytically resummed transverse 
momentum distribution to the ones obtained with a parton shower using the POWHEG 
method in Ref. ^14J and find significant differences for small transverse momenta when 
considering both top- and bottom-mass effects. 

The remainder of the paper is organized as follows: In the next section we will introduce the 
formalism we used for the transverse momentum resummation, indicate the resummtion 
coefficients for Higgs production in gluon fusion and discuss the inclusion of the top- and 
bottom-mass effects in the resummed cross section. We outline our calculation in Sect.[3l 
In the final part of the paper we discuss the mass effects on three quanties: After recalling 
the known results for the total cross section and the LO transverse momentum distribution, 
we present the resummed transverse momentum distribution of the Higgs including the 
full mass dependence at LO+NLL. We conclude in Sect.[5| 



2 Transverse momentum resummation 



2.1 Method 



In this section we sketch the formalism developed in Ref. 29 to resum large logarithmic 



contributions of the transverse momentum (pt) of a colorless particle in the final state. 
We stick to the case of a Higgs boson which is produced in gluon fusionj^ We refer to 
Ref. |29| and references therein for a more detailed description. 



Due to collinear and soft singularities a fixed order transverse momentum distribution 
diverges at small pT- This divergence is evident in the logarithmic pT structure of the 
fixed order (f.o.) cross section 



da 



f.o. 



dp"^ 



TT 



Pt 



In 
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Pt 



Pt 



+ x('-^'^) + 0{pl/ml)]+Oial 



(1) 



The formalism in Ref. |29| is based on the all-order transverse momentum resummation method 



developed in Ref. 31-38 . 
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where rriH and as denote the Higgs mass and the strong couphng constant, respectively. 
The X^i;'*) (a = 1,2) are introduced as the logarithmic coefficients and X^-^'^^ as the 
constant term in p'^ at LO. To obtain a reliable cross section prediction for small transverse 
momenta, these logarithms have to be resummed to all orders in perturbation theory, 
leading to a finite distribution. In the notation of Ref. 39 the hadronic formula for the 
resummed logarithms integrated over the rapidity of the Higgs (y) read^ 



dp, 



a;b={q,q,g} 
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where the Sudakov factor is given by 
S{rn\j, h) = exp < — 



H dq^ 
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The coefficients A, B, C and H are defined as a power series in ag 



X2 tjQ 
Z2'b^ 



(2) 
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Cab{z)=Sab6il-z)+J2 
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(4) 
(5) 
(6) 



A, B and C are process independent, but depend on the initial states of the LO process 
{gg or qq). Accordingly, as soon as these coefficients are determined for one process, they 
are determined for all processes induced by the same initial states. The only dependence 
on the process is embodied in the coefficient once a particular resummation scheme is 
choosen f4l1|^ Instead of H one can also determine the so-called partonic hard-collinear 
function H |39|, introduced in Ref. [29j, which is defined through C and H 



Ugg^abiz) = H / dzi / dz2 5{z - ZiZ2)Cga 

Jo Jo 



(Zl) Cgb{Z2), 



(7) 
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here and in what follows, since it enters at 



Note that we omit the coefficient G introduced in Ref. 
NLO+NNLL which is beyond the accuracy needed in this paper. Furthermore, there is a subtle importance 
of different arguments of Os in the original formula which is not expressed in this formula, since it is not 
essential for the purpose of this paper. 

''Note that only when fixing the resummation scheme the coefficients H, B and C are unambiguously 
defined, since they are connected through so-called resummation-scheme transformations. Fixing H (or 
C) for a single process amounts to fixing the resummation scheme (41]. 
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which again can be written as a power series in as 

oo 
1=1 

Let us make a few more comments about the resummation formula in eq. ([2]): 

• The transverse momentum resummation has to be carried out in the impact pa- 
rameter space (6-space) instead of pT- The integration over b is basically an inverse 
Fourier transform from b- to pT-space. Jo{bpT) is the Oth-order Bessel function and 
the numerical coefficient 60 is the Euler number bo = je = 0.577... . 

• fj'-'^^ denotes the LO cross section of the partonic process gg — )• 

• The integration limits of the parton fractions zi and Z2 are xi = niH / ^/s and 
X2 = mill yfs^ where s denotes the hadronic centre-of-mass energy. 

• The functions fi denote the parton distribution functions (PDFs). 

• For practical reasons (see, e.g., Ref. [29] ) transverse momentum resummation is 
usually performed in Mellin space (iV-space). For this purpose one has to use the 
A''-moment^of all quantities with respect to z = m|^/s, where s denotes the partonic 
centre-of-mass energy. 

• da^'^^ / dp^ in eq. ([2]) fulfills the unitarity constraint, which is imposed in eq. (8) of 
Ref. j29]. More precisely the coefficient % can be choosen in such a way that the 
integral over px of the matched cross section, defined below, reproduces the total 
cross section. This serves as an important consistency check of the matched cross 
section. 



Let us assume that the corresponding resummation coefficients of a process are known. 
Then we have to match the fixed order cross section in eq. ([T]), which is valid for high 
transverse momenta, and the resummed logarithmically enhanced contributions in eq. ([2]), 
which describe the cross section in the small-pT region, to obtain a continuous transverse 
momentum distribution. For this purpose we subtract the logarithmic terms truncated at 
a fixed order 

dpljp 

^ Since the LO of the process gg ^ H does not correspond to the LO of the pt distribution of the Higgs, 
we will refer to ct^"^ as the Born factor in the following to avoid confusion. 
®The A'^-moments of a function g[z) are defined as t/jv = dz z^~^ 9{^)- 



da' 



dpi. 



f n 



(9) 
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from the cross section in eq. ^ at the same order. After that, we add back the logarithms 
resummed to all orders at the corresponding logarithmic accuracy (La.). 

\dpj, J dpj, dp J, 

This features two properties: First, the divergence of the fixed order cross section is 
canceled, leading to a finite result on the right hand side of eq. ( |10[ ). Furthermore, double 
counting of logarithmic contributions present in both da^'°' /dp\ and da^°^ /dp^ is avoided. 

In the formalism of Ref. [29j a new scale is introduced, namely the resummation scale 
Qxcs- For simplicity we refrained from including it into the formulas of this section and set 
it equal to the Higgs mass. This scale may serve as an indicator of the uncertainty induced 
by the truncation of the cross section at some logarithmic accuracy. We will use it together 
with the other unphysical scales, namely the factorization [fj-p) and renormalization scale 
(/i/j), to obtain an error estimate of the cross section. 



da' 



dp% 



1 a 



(10) 



2.2 Gluon fusion at LO+NLL 

In this paper we consider transverse momentum resummation of the gluon fusion process 
at LO+NLL, where the Higgs is produced through both a top- and a bottom-quark loop. 
While higher order corrections are usually evaluated in the heavy-top limit, our calculation 
includes the full dependence on the top {mt) and bottom mass (m^). 



To calculate the resummed cross section in eq. (10) one has to know the resummation 



coefficients described in Sect. 2J^ It is sufficient to determine the coefficient A^'^^ to resum 
all logarithms at leading logarithmic accuracy. For the full NLL resummation A^^\ B^^^ 
and are needed, while at NNLL A^'^\ A^'^\ B^'^^ and Ti^'^^ are required in addition 

For Higgs production in gluon fusion the NLL coefficients ar^ 

A^'^ = Ca, (11) 
BW = -l3, = -(^^CA-lnjy (12) 

^£U = ^(i--)(c^^y + ^-^)' (13) 

^£l,,(z) = ^SU,(^) = -l Pg,iz) = \Cfz, (14) 
= 0, (15) 



^SeeRef. [29j. 

*We set Hf = IJ-R ~ Qrcs = tuh throughout this chapter. 
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where = 5 is the number of hght quark flavors, Ca = 3 is a constant and A denotes 
the finite part of the virtual correction^ which makes T-L^ggi^ggiz) process dependent. 

2.3 Mass effects 

One of the main goals of this paper is the comparison of the resummed cross section with 
exact top and bottom masses to the heavy-top limit. Let us illustrate how the mass effects 



enter the cross section in eq. (10). 



As obvious from eq. ^ and ^ both da^'^^/dp'^ and da^^^^/dp^ are proportional to the 
Born factor a^^\ Thus these contributions calculated in the pure heavy-top limit can 
be reweighted by simply replacing the Born factor calculated with an infinite top mass 
"^hti Born factor with the full top-mass dependence 0"^*^^ or with the full top- and 

bottom-mass dependence (t[^-^. After reweighting, there is no difference between those 
contributions in the heavy-top approach and those including the full mass dependence 
except for the resummation coefficient Ti^^^ of the ^ig-channel, where the mass effects 



enter through the finite part of the virtual corrections A, see eq. (13). 

All other mass effects concern the fixed order contribution da^'°' / dpT of the resummed cross 
section, which enter already at the amplitude level. Nevertheless, in the limit pT ^ 
the fixed order cross section factorizes as well into the Born factor a^'^^ and a process 
independent splitting function. As a consequence, at fixed order the reweighted heavy-top 
limit and the cross section with full mass dependence have the same small-pr limit. This is 
essential and actually a mandatory condition, since the mass effects in da^°^^ /dp^ factorize 
entirely into a^'^\ see above, and it has to subtract the divergence of the fixed order cross 
section for small transverse momenta. 

Since it will be needed throughout this paper, let us define the reweighted cross section as 
follows: 

dax-^Y = dax ■ (16) 

E.g. cjhti-s>t+b = fhtr ('^t+b/'^hti) total cross section in the heavy-top limit reweighted 

by the Born factor including the full top- and bottom-mass dependence. In particular 
daht\-^t+h / dpT is the transverse momentum distribution of cThti->-t+b- When the exact mass 
effects are not known, the reweighted cross section serves as a possible approximation of 
them. While the top-quark effects can usually be well described by the reweighted cross 
section in the heavy-top limit, it is clear that in general this is not true for the bottom- 
quark contribution. We will study this statement quantitatively in Sect.|4j 



'We define the finite part of the virtual according to eq. (38) of Ref. 



42 
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(a) (b) 

Figure 1: (a) LO Feynman diagram of the gluon fusion process and (b) sample 
Feynman diagram of the virtual corrections. The graphical notation for the lines 
is: Thick straight = top or bottom quark; spiraled = gluon; dashed = Higgs 
boson. 



3 Outline of the calculation 



Considering Higgs production in gluon fusion only the diagram shown in Fig.[T](a) has to 
be taken into account at LO. The top quark gives the dominant contribution, since the 
Yukawa coupling is proportional to the corresponding quark mass. While the bottom loop 
still has a considerably large impact on the cross section, all other quarks can be safely 
omitted. 

At NLO the virtual corrections, see Fig.[l](b), and the real emission diagrams in Fig. [2] 
have to be included. To avoid the approximation of an infinitely heavy top quark, all 
contributions need to be calculated by keeping the full top- and bottom-mass dependence. 
The total cross section has already been evaluated with exact top and bottom masses at 



NLO some time ago [19], see also Ref. |4S-47 



In Fig.[T] the Higgs' transverse momentum is always zero, since it is the only particle in 
the final state. Therefore, only the real emission diagrams in Fig.[2j which contribute to 
the total cross section at NLO, enter the transverse momentum distribution at LO for 
Pt > 0. Accordingly, these diagrams contribute to the first term on the right hand side of 



the resummed cross section in eq. ( 10 ) at LO+NLL: 



da 
dp\ 



LO+NLL 



da 



LO 



dpi 



da' 



dpi 



+ 



LO 



da' 



dpi 



(17) 



NLL 



As described in Sect.[2j there are two more contributions to the resummed cross section 
that contain mass effects: The Born factor, which enters both in the second and third 



term on the right hand side of eq. (17); and the finite part of the virtual corrections. 



which enters in the coefficient H^^^ of the gg-channel, see eq. (13). These contributions 



with full mass dependence, the Born factor, the virtual and the real matrix elements, were 
obtained from the authors of Ref. [48]. We implemented them into a numerical program 
that calculates the resummed transverse momentum distribution at LO+NLL analogous 
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(a) 



(b) 



(c) 



Figure 2: Sample Feynman diagrams of the real corrections for the individual 
channels (a) gg (b) gq and (c) qq, contributing to the process pp ^ H at NLO 
QCD. Notation as in Fig. [T] and thin straight = light quark q £ {u, d, c, s, b}. 



to HqT [28-30 . The difference is that HqT covers only the heavy-top limit, but up to 



NLO+NNLL. The third term on the right hand side of eq. (17), which resums the large 



logarithmic contributions, was obtained using a modified version of HqT. 

The resummation can be done separately for the individual sub-channels gg (Fig.[2](a)), 
gq (Fig.[2](b)) and qq (Fig.[2](c)). The g^r-channel contains double and single logarithmic 
contributions, while the gg-channel only embodies a single logarithmic divergence at small 
transverse momenta. The gq-channel on the other hands remains finite for small px on 
its own, since it contains no logarithmically enhanced terms at LO. Consequently no 
resummation has to be performed in this channel. 



Since the resummation procedure of Ref. [29] fulfills the unitarity constraint^ we were 
able to use this as an important cross-check of our calculation. For the integral of the 
resummed cross section over all pT, we found agreement at the sub-percentage level with 
the total cross section [l9], separately for the individual sub-channels. We also checked 
that this integral is resummation scale independent. Furthermore, we compared the fixed 
order transverse momentum distribution to Ref. [48] and checked all contributions in the 
heavy-top limit against HqT. In both cases we found perfect agreement. 



4 Results 

4.1 Total cross section 

Before we discuss the transverse momentum distribution with exact top- and bottom-mass 
dependence, let us first summarize the quark-mass effects in case of the total inclusive cross 
section. 

In the heavy-top limit, the total inclusive cross section has been calculated at NNLO [4]-[6], 
i°See Sect.ET 
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while the full quark-mass dependence is only known through NLO |19j. One observes that 
the difference between the reweighted heavy-top limit and the exact top-mass dependence 
at NLO is at the sub-percentage level for a Higgs mass lower than twice the top mass. 
In particular this means that the corresponding ratios of the NLO and LO cross sections 
Khti = ^hti°/'^hti ~ '^^^^'^ / '^t'^ (iT- factors) agree within this accuracy. A few years 

ago, finite top-mass effects have been investigated also at NNLO (8- 10 , showing that these 
effects remain below 1%. 

When the bottom-quark contribution is included, the total inclusive NLO cross section 
is reduced by about 7% for a Higgs mass of 125 GeV at the LHC with 8 TeV machine 
energy. At LO, however, the cross section is lowered by more than 10%. Evidently, the LO 
bottom-quark contribution is negative while the effect on the NLO correction is positive. 
Hence the i^-factors ii'hti and Kt differ from -fCt+b a-t the order of 3.5%. 

If the bottom contribution was not known at NLO, one could try to approximate it through 
reweighting of cr^^^^ by Kt or K^ti- A comparison with the exact NLO cross section cr^^^ 
shows that it is better to not reweight the NLO cross section and therefore to omit the 
NLO bottom-quark contribution completely, but take into account only the exact LO 
dependence. 

A similar conclusion might be true for differential quantities as well. Therefore we present 
the full quark-mass dependence of the resummed transverse momentum distribution in 
this paper and investigate the mass effects. 



4.2 Transverse momentum distribution 
4.2.1 Preliminary remarks 

After some considerations concerning the pT distribution at LO, we will show the resummed 
cross section with respect to the transverse momentum of the Higgs at LO-I-NLL including 
the full top- and bottom-mass dependence. We present results for the LHC at 8 TeV and a 



Higgs mass of 125 GeV. We use the MSTW2008 NLO PDF sets [49] and the corresponding 
input of the strong coupling constant as{mz) = 0.12018, where mz denotes the Z-Boson 
mass. Our central scale choice is /io = nT-Hf^ for the renormalization, factorization and 
the resummation scale. We insert an on-shell top and bottom mass of rrif = 172 GeV and 
nib = 4.9 GeV, respectively. 



4.2.2 LO considerations 

Before analyzing the quark-mass effects on the resummed transverse momentum distribu- 
tion at LO+NLL, let us recall the situation for the pure LO px distribution. It is shown 
in Fig. ^ for the heavy-top limit da^^^/dpT (black, dotted), exact top-mass dependence 
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Figure 3: LO transverse momentum distribution of the Higgs in the heavy- 
top limit (black, dotted), including the exact top mass (blue, dashed) and with 
exact top- and bottom-mass dependence (red, solid). The uncertainty bands are 
obtained by scale variation. In particular we varied one scale within [0.5 fiQ, 2 ^o]> 
while fixing the other, and vice versa. The maximum and minimum values of the 
cross section in this procedure are the uncertainties to the central value. In this 
logarithmic plot the curve of dat/dpT is almost indistinguishable from dat^^/dpT- 



dal''^/dpT (blue, dashed) and including both top- and bottom-quark masses da^^^/dpx 
(red, solid). As expected, the fixed order curves diverge in the limit pT — >■ 0. Further- 
more, all three curves are quite close to each other for small transverse momenta, while for 
Pt > 200 GeV a gap between the heavy-top limit and the other two curves emerges. Since 
the top contribution is dominant, daf^ /dpx and da]^^^/dpT stay close to each other as well 
for high transverse momenta. The uncertainty bands are obtained through independent 
variation of and within [0.5/iOi 2 

Considering the relative contribution of the mass effects normalized to the heavy-top 
limit, see Fig.|4j we can visualize discrepancies also for small transverse momenta. For 
comparison we show three different curves for the heavy-top limit (black, dotted): The 
pure heavy-top limit is the normalization in this plot and hence located at one. The other 
two lines denote the reweighted cross sections dal^^^t/'^PT (upper line) and dal^2^^.^y^/dpT 
(lower line), respectively. 

Let us first compare the cross section with full top mass dependence (blue, dashed) to 
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Figure 4: Curves of Fig. [3] normalized to the heavy-top hmit curve. The upper 
and lower black dotted curves denote the reweighted cross sections da^2-^t/^PT 
and dal^2^^+h/^PTj respectively. 



the reweighted cross section da^^^^-^^/ dpT- We find that, similar to what was found in 



Ref. 11 , the reweighted heavy-top limit works for the top loop to better than 4% as 
long as pt < 150 GeV. For large transverse momenta the heavy-top approximation is 
not valid any more, because the discrepancy is growing rapidly for px > 150 GeV. The 
kink at pT ~ 150 GeV occurs in the qq channel and is due to the kinematical cut at 



/p^ + mjj + Pt- Because of the factorization of the cross section in the limit 
Pt — )• 0, see Sect. 2.3, da\'^ /dpT and da^^^^^/dpT in Fig. |4] become identical in this limit. 



The same is true for the cross section including top- and bottom-mass dependence (red, 
solid) in Fig.|4j da\2\^/dpT converges to da^^Y_^^j^^^/ dpT in the limit pT — )■ 0. Let us compare 
the cross section with exact top and bottom masses to the pure heavy-top limit (black, 
dotted) located at one. The discrepancy ranges between —5% and +7% for p^ < 200 GeV. 
The difference is even higher (up to 13%) when we compare da^^^^/dpx to the reweighted 
cross section dal^^^^_^_^/dpT- Evidently, as soon as the bottom loop is considered, omitting 
the mass effects works better than approximating them by the reweighted cross section 
^^hti^t+h/ '^Pt ■ This is similar to the total cross section, where the NLO i^-factors -fChti 
and Kt should not be applied to the bottom-quark contribution. In the following we will 
investigate if a similar conclusion can be drawn for the resummed transverse momentum 
distribution. 
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Figure 5: Same as Fig.[3j but for the resummed cross sections and splitted into 
two ranges, from 0-100 GeV and 100-500 GeV. 



4.2.3 Resummed cross section at LO-|-NLL 

Fig. [5] shows the matched cross sections of the fixed order prediction and resummed large 



logarithmic contributions at LO-I-NLL, defined in eq. (17). As before we plot three different 
curves, the pure heavy-top approximation (black, dotted), the cross section including full 
top-mass dependence (blue, dashed) and the one with exact top and bottom quark masses 
(red, solid). Unlike the LO distribution, the resummed cross sections are finite at small 
transverse momenta, leading to a reliable prediction also for pT ^ 30 GeV. 

The uncertainties of the cross sections emerge from the truncation of perturbative series 
with respect to both the strong coupling constant and the logarithmic accuracy. They are 
estimated by the variation of the renormalization and factorization scale and the resum- 
mation scale, respectively. The error bands are obtained through independent variation 
of iiR, lip and Qres within [0.5 ^uq, 2 ^uq]- We neglect the uncertainty emerging from the 
PDFs and a^. It was already found for the heavy-top limit, see Ref. (29j Fig. 8, that for 
Qres ^ w-J? the matched cross section can become negative for high transverse momenta. 
Accordingly, the uncertainty bands of all cross sections in Fig. [5] become negative around 
Pt ~ 150 GeV because of the resummation scale variation. For pT ^ 100 GeV the relative 
scale uncertainty ranges between 10% and 50%, while in average they are of the order of 
25%. For higher transverse momenta the uncertainty gets huge, due to the fact that the 
error bands become negative. 

Let us investigate the mass effects on the resummed cross section in the low px range 
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Figure 6: Same as Fig. |4] but for the resummed cross sections. 



in more detail. For this purpose the relative contributions normalized to the heavy-top 
limit are plotted in Fig.|6j Just as in Fig. [4] there are three curves for the heavy-top limit 
(black, dotted). The one in the middle denotes the pure heavy-top limit, while the other 
two denote the reweighted cross sections. The upper line represents da^^^^^^/dpr and 
the lower one dal;^:^^^^ / dpx ■ 

Let us first compare the cross section with full top-mass dependence (blue, dashed) to the 
reweighted cross section da^^^^^^ / dpx ■ We find that the reweighted heavy-top limit is 
working to better than 4.5% as long as px < 150 GeV, while, as expected, the discrepancy 
grows as the transverse momentum of the Higgs increases. In the region pT < 50 GeV, 
where resummation effects are important, the top-mass effects remain even below 0.5%. 
In the small-pT limit the resummed cross section including the exact top-quark mass is 
not identical to da^^^^^^^ / dpr in Fig.|6]as it was for the LO distribution. This is caused 

by the fact that the mass effects do not completely factorize into the Born factor see 
Sect.EH 

The same is true for the resummed cross section with full top- and bottom-mass de- 
pendence (red, solid) in Fig.[6| where the gap to da^^^^^^/dpx in the limit px ^ is 
even larger. Let us investigate the mass effects when both top and bottom quark are 
included and determine the quality of the reweighted cross section da^^^^^^^ / dpx to ap- 
proximate them. Therefore, we compare the red, solid curve {da^^^^^^ / dpx) to the lower 
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Figure 7: Resummed cross sections normalized to the one including the full 
top- mass effects. 



black, dotted line {da^^^l^^^ / dpx) and find that the discrepancy reaches up to 14% for 
Pt < 200 GeV. In fact omitting the mass effects completely works better: The difference 
between da^^-^^^^ / dpT and the pure heavy-top limit ranges only from —2.5% to 7% for 
Pt < 200 GeV. Note also that omitting the 6-loop contribution and account only for the 
top-mass effects through dal;^yij^^^ / dpx is closer to da^;^^^^^ / dpr than rfcr^°-it+b/^^'7^' 
since in this case the difference reaches only up to 9%. 

Regarding this, it is interesting to examine the bottom-mass effects assuming the exact 
top-mass dependence is known. In Fig. [7] the ratio of the resummed cross section including 
top- and bottom-mass dependence and the one including only top-mass effects is shown 
(red, solid). Omitting the 6- loop effects amounts to about ±8% when the top mass is 
taken into account exactly. 



In an earlier study [14| the mass effects on the resummed transverse momentum distri- 
bution were calculated using the POWHEG method in combination with the PYTHIA 
parton shower |50j. For comparison we reproduced the lower right plot in Fig. 3 of 
Ref. [14| in the case of the analytically resummed cross section at LO+NLL, see Fig.jsj 
It shows dat/dpT normalized to dahti~>-t / dpr (blue, dashed) and d(Tx_j^\^/ dpx normalized 
to dcThti-s.t+b/'^PT (red, solid). Although the parameter choices of Ref. [l4| are slightly 
different, for px ^ 200 GeV the top-mass effects on the analytically resummed cross sec- 
tion are in good agreement with the ones of the resummed cross section obtained with 
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Figure 8: Analogous to lower right plot in Fig. 3 of Ref. [M]. Blue, dashed 
curve: Ratio of dat/dpT and daiiti^t/dpr', red, solid curve: Ratio of dat+h/dpr 
and dam-^t+h/dpT- 



POWHEG+PYTHIA (blue, dashed in Fig.|8j black, dashed in Fig. 3 of Ref. [W]). Considering 
an exact top- and bottom-mass dependence (red, solid curve in our plot; blue, solid curve 
in Ref. [l4]), on the other hand, the mass effects on the two approaches appear to be 
considerably different for px < 50 GeV, i.e. in the region where resummation becomes im- 
portant. However, both approaches are theoretically well defined and the numerical results 
are consistent within the respective resummation formalism. The discrepancy might be 
caused by the different treatment of the next-to-leading logarithms in the two approaches. 
Furthermore, it needs to be clarified whether the discrepancy arises from the normaliza- 
tion factor, or whether it is a genuine effect in the cross section with full mass dependence. 
Currently, it has to be considered as a measure of the theory uncertainty at small pT- 
Clearly, the source of the difference deserves further investigation!^ 

In summary we find that in all cases studied in this paper it is not a good approximation 
to account for 6-loop effects using the reweighted cross section in the heavy-top limit. We 
conclude that bottom-mass effects should be included only up to the order where their 
calculation is feasible and should be omitted otherwise. 



According to Ref. 51 in the MC@NLO approach [52 the shape of the curve including top- and 
bottom-mass dependence is much more similar to ours (red, sohd curve in Fig.js]). 
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Along these lines we give a recommendation for the evaluation of the best prediction 
regarding the resummed cross section of the gluon fusion process. As stated above, trans- 
verse momentum resummation of the Higgs boson is known in the heavy-top limit at 



NLO+NNLL 128-30 , while in this paper we presented the full top- and bottom-mass de- 



pendence at LO+NLL. Both cross sections should be combined by 



est 



J^NLO+NNLL 



dpi 



dpi 



J^LO+NLL j_LO+NLL 
"^htl _^ "^^t+b 



dpi 



dpi 



(18) 



The first two terms on the right hand side of eq. (18) can be calculated with the program 
HqT |28| - (30] . As well for the second but especially for the third term numbers can be 
obtained from the authors of this paper upon request. 



5 Conclusions 



We presented the transverse momentum resummation of the Higgs boson in gluon fusion 
including the full top- and bottom-mass dependence at LO+NLL. We found that top-mass 
effects can be approximated in the heavy-top limit to better than 0.5% for pj- < 50GeV 
and 4.5% for pT < 150 GeV, respectively. Furthermore, we analyzed the contribution of 
the 6-loop in the resummed cross section. It amounts to about 10%. In comparison with an 
earlier study, we found that for small transverse momenta the influence of the bottom mass 
on the analytically resummed transverse momentum distribution is considerably different 
from the one obtained using the POWHEG method in combination with the PYTHIA parton 
shower, while the top-mass effects on both approaches are in good agreement. We also 
showed that the bottom-mass effects are not well approximated by the reweighted cross 
section in the heavy-top limit. As a consequence we recommended a best prediction for 
the resummed cross section that combines two calculations done with current technology: 
The LO+NLL contribution should be subtracted from the NLO+NNLL prediction in the 
heavy-top limit and added back including the full top- and bottom-mass dependence. 
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